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ABSTRACT 

We investigate the ^-process in an energetic Type Ic supernova (SN Ic) and resultant productions of 
the light elements including boron and its stable isotopes. SN Ic is a very uni que boron source because 
it can produce boron through not only spallation reactions as discussed in INakamura 8z Shigevamal 
(2004) but also the ^-process. The ^-process is considered to occur in core-collapse supernovae and 
previous studies were limited to Type II supernovae (SNe II). Although the progenitor star of an SN 
Ic does not posses a He envelope so that 7 Li production via the v- process is unlikely, n B can be 
produced in the C-rich layers. We demonstrate a hydrodynamic simulation of SN Ic explosion and 
estimate the amounts of the light elements produced via the z/-process for the first time, and also the 
subsequent spallation reactions between the outermost layers of compact SN Ic progenitor and the 
ambient medium. We find that the ^-process in the current SN Ic model produces a significant amount 
of n B, which is diluted by 10 B from spallation reactions to get closer to B isotopic ratios observed 
in meteorites. We also confirm that high-temperature /i and r neutrinos and their anti-neutrinos, 
reasonably suggested from the compact structure of SN Ic progenitors, enhance the light element 
production through the neutral-current reactions, which may imply an important role of SNe Ic in 
the Galactic chemical evolution. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: general — stars: 
abundances 



1. INTRODUCTION 

Stars with mass more than about 10 Mq are known to 
culminate their evolution as supernova explosion. Most 
of the gravitational energy released from formation of 
proto-neutron star is carried away by neutrinos and total 
number of emitted neutrinos amounts to more than 10 58 . 
Although cross sections of neutrino-nucleus reactions are 
very small, a fraction of such a huge number of neutrinos, 
in particular \x- and r- neutrinos with higher tempera- 
tures than electron neutrinos, can scatter off nuclei to 
the excited states, induce ejection of neutron, proton, 
or a-particle, and leave observational signature of the 
^-process nucleosynt hesis in the yields of some species 
(jWooslevet al.l[T990T) . 

Massive stars are one of the astrophysical sites for the 
production of light elements via the ^-process, mainly 
contributing to the production of 7 Li and 11 B (e.g. 
Wooslev et all [19901 iHeger et al.l [20051 lYoshida et at] 
20041 ). 7 Li is mainly produced in He-rich envelopes and 
n B in inner C-rich shells. The stellar materials includ- 
ing newly-synthesized light elements are accelerated by 
strong shock wave, a part of which is destroyed by shock 
heating, and finally ejected into the interstellar matter 
to make a partial contribution to the next generations of 
star formation. 

Light element isotopes are also known to be pro- 
duced from some other processes; 7 Li from big bang 
nucleosynthesi s, 7 Li and n B from AGB stars and no- 
vae (|Cameronl 119551 ). and all stable isotopes of Li, Be, 
and B (LiB eB) from Galactic cosmic-ray (CR) nucleosyn- 
thesis (e.g. iMeneguzzi et al.lll97lUSuzuki fc Inouell2002l: 
IRollinde et al.l 120081) . 



Despite a great deal of efforts dedicated to investigate 
these processes, the origins of LiBeB are not fully un- 
derstood and there remain some observational features 
which conflict with theoretical predictions. For instance, 
theoretical predictions for LiBeB production from Galac- 
tic CRs alone indicate quadratic relation between BeB 
abundances and metallicity, while observations clearly 
show linear relation . CNO CRs from superbubbles 
(Higd on et al.l Il998h . supernovae or Wolf-Rayet stars 
(|Ramatv et al.l 119971: lYc-shii et all I1997T). and Type Ic 



super novae (jFields et all 12002b INakamura fc Shigevamal 
I2004D have been suggested to account for the observed 
linearity. However, the CR spallations cannot be the 
only source of boron isotopes because the CR spallations 
do not reproduce the high 11 B-to- 10 B abundance ratio as 
observed in meteorites, and another source of n B such 
as supernova (/-process is necessary. 

It is required to accurately and quantitatively evaluate 
each contribution of the LiBeB production process in or- 
der to solve these problems. In this Letter we report that 
energetic core-collapse Type Ic SNe (SNe Ic) could be a 
viable astrophysical site where light elements including 
boron isotopes are produced. Our nucleosynthesis model 
takes account of the ^-process in SNe Ic for the first time 
and also spallation reactions of accelerated SN ejecta in- 
teracting with interstellar/circumstellar matter, which 
do not conflict with the observed linear relation between 
BeB abundances and the metallicity in metal-deficient 
halo stars. The progenitor of an SN Ic is a C/O star and 
its H and He envelopes have been stripped during the 
stellar evolution. Although the explosion mechanism of 
SNe Ic has not been clarified, neutrinos should be one of 
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the main carriers of the gravitational energy release from 
the collapsing core. The robustness of neutrino emission 
from a collapsing proto-neutron star which was evolved 
from a ~ 40M rn progenitor star was studie d in numeri- 
cal simulation (jSumiyoshi et al.l l2006. 2003). Even in the 
case that the central core of such a star eventually turns 
into a black hole, a temporally formed proto-neutron 
star is shown to emit a huge amount of neutrinos be- 
fore collapsing to the black hole. Even after the black 
hole formatio n, neutrinos are be ing emitted from accre- 
tion disk (e.g. iKohri et" al. 2005; Surman & McLaughlin 
2005). Therefore, a huge amount of neutrinos are emitted 
from the central region of an SN Ic so that the (/-process 
can operate in the exploding supernova material. In §2, 
we present our model of supernova explosion fi )2.1l) . the 
^-process nucleosynthesis ( §2.2|> . and the spallation re- 
actions ( §2.3|) . We discuss our calculated results in §3 
where we clarify the important role of SNe Ic in the light 
element synthesis. 

2. CALCULATIONS AND RESULTS 

2.1. Supernova Explosion 

We consider a very energetic explosion of a 15 M C/O 
star with the explosion ene rgy E ex = 3 x 10 52 er gs cor- 
responding to SN 1998bw (jNakamura et all 120011) . The 
explosion energy is released at the center of the progen- 
itor star as thermal energy. Resulting shock wave accel- 
erates the stellar materials and explodes the progenitor 
as a supernova. Time evolutions of physical quantities 
in the progenitor are calculated with ID hydrodynamic 
code which takes account of the effects of special relativ- 
ity. We solve the special relativistic hydrodynamic equa- 
tions in Lagrangian coordinates with an ideal equation 
of state involving gas and radiation pressure. Adiabatic 
indices are treated as functions of pressure and gas den- 
sity. Details on the numerical method are described in 
iNakamura fc Shigevamal (|2004D . 
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Fig. 1. — The mass fraction distributions of light elements 
produced via the i/-process as functions of the mass coordinate. 
Shown are the cases of the SN 1998bw model (M co = 15M , 
E cx = 3 X 10 52 ergs) with T„ e = 3.2 MeV, T,7 C = 5 MeV, and 
T„ MT = 6 MeV (top) and 8 MeV (bottom) at 75 seconds after the 
energy release. 



2.2. The Neutrino- Process 

Light element synthesis in the Type Ic supernova is 
calculated as a post process. We use a nuclear reaction 
network consisting of 29 1 species of nuclei a nd taking 
account of the ^-process (jYoshida et al.ll20C)8J) . 

Neutri nos are emitted from a collapsing proto-neutron 
star (e.g. lSumiyoshi et al.ll200 6) and /or the innermost re- 
gion j ust above a black hole (e.g. Surman & McLaughlin 
2005). Neutrino properties, particularly in SNe Ic, are 
still un certain. We therefo re use a supernova-neutrino 
model ([Yoshida et al.l [2008h that the neutrino luminos- 
ity decreases exponentially with a time scale of 3 s and 
that the neutrino temperature of each species does not 
change with time. The total energy carried out by neu- 
trinos is assumed to be 3 x 10 53 ergs. 

We consider two different models for neutrino temper- 
atures. For both models the temperatures of v e and v e 
are set to be T Ve =3.2 MeV and T Pe = 5 MeV. In the 
first model, the temperature of the other flavors, ^. T 
and Vfi,Ti takes an identical value of T v T =6 MeV. This 
temperature model is referred to as the "standard" T v 
model. The light element synthesis in an ~ 2OM Type 
II supernova with the standard temperature model well 
reproduces the supernova contribution of t he 11 B produc- 
tion during Galactic chemical evolution (Yoshida et al. 



120051 [20081 ). The third peak of r- process elements is also 
reproduced well in neutri no-driven win d models using 
these temperatures (| Yoshida et a l. 2004). In the second 
model, the temperature of z/ MjT and v^^ is set to be T„ 
= 8 MeV. This temperature model is referred to as the 
"high" T v T model. A newly forming proto-neutron star 
in an SN Ic should be more compact than that of an SN 
II. Therefore, the neutrino temperature of the collaps- 
ing proto-neutron sta r is larger than that of SN II (e.g. 
iSumivoshi et~all[2006l ). 

Calculated mass fraction distributions of n B, n C, 10 B, 
7 Li, and 7 Be in the SN Ic model are shown in Figure 
1. The yield of each species strongly depends on the 
branching ratios of the decay channel of 12 C-neutrino 
reactions. The light elements are mainly produced in the 
O/Ne layer (M r > 7AM ). The mass fractions in the 
O/Si layer (4.6M Q < M r < 7AM Q ) are smaller than in 
the outer layer. The Si/S layer (3.7M Q < M r < 4.6M ) 
where carbon completely burned out presents no light- 
element production. The mass fractions of some light 
elements are also large in the innermost region where en- 
rich freeze out is achieved (M r < 2.6M Q ) and complete 
Si-burning occurs (2.6M < M r < 3.5M ). There is no 
qualitative difference in the mass fraction distributions 
between the standard and high T u T models except that 
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TABLE 1 
Light Element Yields from SN Ic Model 



Species 


Standard T„ T 


High T„ MiT 


Spallation 




(M ) 


(Me) 


(M@) 


6 Li 


1.67 x 1CT 11 


5.91 x 10~ n 


2.38 x 10~ 7 


7 Li 


7.41 x 10~ 9 


2.50 x 10~ 8 


3.31 x 10" 7 


9 Be 


4.49 x 10" 11 


1.08 x lO" 10 


9.99 x 10" 8 


io B 


1.29 x 1CT 9 


2.78 x 10~ 9 


4.38 x 10~ 7 


n B 


2.69 x 10~ 7 


5.46 x 10~ 7 


1.34 x 10~ 6 



the mass fraction of each species in the high T v model 
at a given mass coordinate is roughly twice as large as 
the one in the standard T v r model. 

Several light elements produced from 12 C- neutrino re- 
actions are exposed to explosive nucleosynthesis when 
the shock arrives. n B, 11 C, and 10 B are destroyed sub- 
sequently by collisions with protons and a-particles. 7 Li 
and 7 Be are photodisintcgrated to 3 H and 3 He, respec- 
tively. After the shock passage, the exploding materials 
are still being irradiated by neutrinos so that light ele- 
ments are produced through the ^-process again. This is 
the reason why larger mass fractions of n B, 'Li, and 7 Be 
in the outermost region (M r > 13. 5M©) survive, being 
avoided from explosive nucleosynthesis. 

In the innermost region, light elements are produced 
after the termination of the nuclear statistical equilib- 
rium. The main product in a-rich freeze out is He. 
About 20 % of 4 He by mass fraction is also produced 
in complete Si-burning. The ^-process of He produces 
3 H and 3 He in cooling materials, followed by a-capture 
reactions to produce 7 Li and 7 Be. Furthermore, n B and 
n C are produced by a-captures, and 10 B by 7 Be(a,p) 10 B 
reaction. 

Table [1] summarizes the total yields of the light ele- 
ments produced through the i/-process. In both models 
of T v significant amount of n B of order 10 -7 Mq is 
synthesized. It should be noted that the amounts of n B 
include n C which decays to n B in 20 minutes. The 
yield of 7 Li is of the order of 10" 9 - 10- 8 M Q in the SN 
Ic model. Thi s yield is much smal ler than that produced 
in SNe II (e.g. lYoshida et al.lf2008h . Most of 7 Li in SNe II 
is produced in the He-rich layer. In contrast, almost all H 
and He layers of SN Ic progenitors have been stripped via 
stellar wind and/or binary effect before explosion. The 
yields of 6 Li, 9 Be, and 10 B are of the order of or below 
1O _9 M0 due to sma ller branching ratios of neutrino- 12 C 
reactions (jYoshida et al. 2 0081 ). 

2.3. Spallation Reactions 

Supernova explosions accelerate and expel the stellar 
materials into the interstellar space, which pollutes the 
universe with newly synthesized elements. The progen- 
itors of SNe Ic are so compact that a small fraction of 
ejecta can be accelerated nearly to the speed of light. 
Thus accelerated surface layers of SNe Ic, which are 
composed of C and O, interact with interstellar matter 
(Nakamura & Shigcvama 2004) or circumstellar matter 
( Nakamura et al.l 120 06?) , and produce the light element 
isotopes via spallation reactions of CNO with protons or 
a-particles. 

Her e we use the m odel for SN 1998bw constructed by 
iNakamura fc Shigevamal IpOOl . About 0.3% (0.04M Q ) 



of the ejecta attain enough energy (> 10 MeV/A) to 
undergo spallations (see eq. [2]). We solve the transfer 
equation for each element i expressed as 



dFi(e,t) _ d[oj t {e)Fi{e,t)] F t (e,t) 



Of 



(k 



A 



f>Vi(e), (1) 



where A is the loss length in g cm' 2 , p denotes the mass 
density of the ISM, Vi(e) the velocity of the element i 
with an energy per nuclcon of e. The initial condition 
for the mass of element i with an energy per nucleon e 
at time t — 0, Fi(e, t — 0), is derived from the numerical 
calculations of explosions or an empirical formula 



^P± = 1.9x10- 



ff cx /10 51 ergs 

M ej /1 M 



3.4 

VIOMeV/ 



-3.6 



(2) 



where M e j is the mass of the ejecta. The yield of a light 
element I via the i + j —> I + • • • reaction is estimated 
from 

dN t f i , ,Fi(e,t) 



dt 



»w- 



AiUip 



■Vi(e)de, 



(3) 



where A) is the number of the produced light element I , 
rij is the number density of an element j in the interstel- 
lar matt er, <j\ 7 - the cros s section of i+j —>/ + ••• reaction 
given bv lRead k. Violal (1984), and A t the mass number 
of the element i. The interstellar matter is assumed to 
consist of neutral H and He with the number densities of 
Tin = 1 cm~ 3 and tiho = 0.1 cm~ 3 . 

Calculated results show that the mass of n B produced 
via spallation reactions is 1.3 x 10~ 6 A/ Q , which is larger 
than that synthesized via the ^-process even in the high 
T v i T model (Table [T]) . The resultant isotopic ratio of 
11 B/ 10 B (~ 3) from spallation reactions in this model is 
predominantly determined by the ratio of cross sections 
of the reaction p,a + O — > 10 B to that of p,a + O — s- u B. 
The yields of the other light elements are listed in Table 

m 

3. CONCLUSIONS AND DISCUSSION 

INakamura fc Shigevamal (|2004h concluded that SNe Ic 
may not play an important role in B isotope productions 
because of the low isotope ratios 11 B/ 10 B ~ 2.8 com- 
pared with observations in meteorites (4.05 ± 0.05). In 
this Letter we estimate the production of light elements 
including boron isotopes via the ^-process in an SN Ic 
and find that SNe Ic can be viable n B producers as well 
as SNe II. The resulting number abundance ratio of B 
isotopes from both the ^-process and spallations turns 
out to be n B/ ln B = 3.7 (for standard T^J -4.3 (for 
high T v ). It is reasonable to assume that the radii 
of neutrino spheres in SNe Ic are smaller than those in 
SNe II because of the compactness of SN Ic progenitors, 
leading to the higher temperatures of emitted neutrinos 
and anti-neutrinos. Higher T„ T results in higher LiBeB 
yields through neutral current reactions in the ^-process, 
which raises the significance of potential roles of SNe Ic 
in the light element production. The 11 B/ 10 B ratios in- 
ferred from current SN Ic models seem to agree well with 
solar values. However, the frequency of SNe Ic, in par- 
ticular with high explosion energy as assumed here, is 
quite low and borons in meteorites should be dominated 
by those originated from SNe II and Galactic CRs. In 
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fact, according to iCappellaro et all (|1999T) the estimated 
frequency ratio of SNe Ib/c to SNe II in local spiral galax- 
ies falls in the range of 0.41 to 0.058 and the neutrino- 
processed 11 B yield from each SN Ic is at most com- 
parable to the contribution from SN II (averaged value 
per su pernova ~ 8.3 x 10~ 7 M from iRauscher et al.l 
(2002); see also the case of SN 1987A mo del with simi- 
lar parameter set in lYoshida et al.1 ((2008)). These facts 
lead to a conclusion that SNe Ic do not dominate the 
bulk of present n B in the Galaxy. Contribution from 
SNe Ic to the light element production in metal-deficient 
stars might be outstanding because an SN Ic progenitor 
is surrounded by its wind material and the light elements 
produced in the explosion are likely to be inherited di- 
rectly by next generations of stars, which would show 
high BeB abundances overlying small amounts of BeB 
from the other processes. An extraord inary Be-rich halo 
star HD106038 (|Smilianic et al.ll2008ft could be a candi- 
date star exhibiting such a process. 

All previous Galactic chemical evolution (GCE) mod- 
els included four different source components of the 
light element productions from primordial nucleosynthe- 
sis, AGB star and nova nucleosynthesis, SNe II ^-process 
nucle osynthesis, and G alactic CR spallations and fu- 
sion (jRyan et al.ll200lL and references therein). Based 
on the GCE s tudies especially of the B isotope ratio, 
lYoshida et al.l (J2005I ) discussed that the energy spectra 
of supernova neutrinos can be constrained through the 
avoidance of an overprod uction of the n B abundance 
(|Wooslev fc Weaver] 119950 . Our present result in this 
Letter indicates that the GCE model should include the 
fifth astrophysical component of SNe Ic where both re- 
process and CNO CR spallation contribute remarkably 
to the production of 7 Li and n B. The GCE studies in- 
cluding both contributions from SNe II and SNe Ic are 
highly desirable. 

Neutrino oscillations change neutrino flavors in stellar 
materials. It has been shown that light element yields in 
SNe II increase when the Mikheyev-Smi rnov-Wol fstcin 
(MSW) effect is taken into account (Yoshida ct al. 
2006a. b, 2008). The contribution from charged-current 
^-process reactions becomes more effective in the He 
layer because the spectra of v e and v e shift to higher 
energy owing to the MSW effect. In the case of SNe 
Ic, however, the increase in the light element yields by 
the MSW effect is expected to be small for the lack of 
He layer. This relates to the locations of two resonances 
of the MSW effect. Large flavor exchange occurs be- 
tween v e (v e ) and z/ PiT (Pfj,, T ) at the higher density reso- 
nance in the case of normal (inverted) mass hierarchy if 



sin 2 26>i3 > 10~ 3 . The density of the higher resonance for 
a neutrino with energy > 20 MeV is smaller than 1,600 
g cm -3 , corresponding to M r > 13.5M Q in the mass co- 
ordinate. Therefore, the ^-process only in the outermost 
region in our SN Ic model will be affected slightly by the 
MSW effect. 

Recently, neutrino oscillations with ne utrino self- 
interactions have been investigated (e.g. iDuan et"aT1 
[20061: iFogliet al.l l2007f ) . The flavor change by this ef- 
fect occurs in deep region of SN ejecta (~ 10 7 cm) and 
the abundances of light element are expected to increase 
in entire region of the SN Ic. This is because high- 
energy part of the neutrino spectra plays major role in 
the ^-process nucleosynthesis where the neutrino self- 
interactions increase the contribution of charge current 
interaction in addition to the invariant neutral current 
interaction in the standard temperature hierarchy of su- 
pernova neutrinos which we set here T Ve < T Pe < T„ T . 
Details of the effects of neutrino self-interactions which 
depend on the initial spectra, luminosity, density profile, 
and mixing angles of neutrinos will be discussed else- 
where. 

The combined mechanism of light element production 
described here, i.e. the ^-process and spallations, may 
also apply to the case of Type lb supernovae where pro- 
genitor stars without H-rich envelope explode. The out- 
ermost layers of SNe lb is composed of He so that the 
a-fusion reaction to pro duce lithium isotopes be comes 
effective. Furthermore, iMevnet fc Maedeii (|2002f ) sug- 
gested that the He layers of very metal-poor stars con- 
tain primary nitrogen which was produced during the 
He-burning phase by the CNO cycle in the H-burning 
shell stimulated by rotationally-induced diffusion of car- 
bon into the H-burning shell. This strong thermonuclear 
reaction may enhance the light element production via 
spallations because the compact progenitor resulted from 
strong mass loss leads to more effective acceleration of 
ejecta by shock wave and the cross sections of spallation 
reactions involving nitrogen show high peak values at 
relatively low energy thresholds. Observationally some 
low-metallicity halo stars were found to have high abun- 
dances of 6 Li or 9 Be beyond the theoretical predictions, 
which engages our interest in SNe lb even if they might 
not accelerate their envelopes effectively because their 
progenitors are not so compact as those of SNe Ic. A vari- 
ety of progenitor models and parameters of SN explosion 
will be considered in the forthcoming paper, including 
their contributions to the Galactic chemical evolution. 
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